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ABSTRACT 

We discuss the systematic uncertainties inherent to analyses of observed (broad-band) 
Spectral Energy Distributions (SEDs) of star clusters with evolutionary synthesis 
models. We investigate the effects caused by restricting oneself to a limited number 
of available passbands, choices of various passband combinations, finite observational 
errors, non-continuous model input parameter values, and restrictions in parameter 
space allowed during analysis. Starting from a complete set of UBVRIJH passbands 
(respectively their Hubble Space Teiescope/WFPC2 equivalents) we investigate to 
which extent clusters with different combinations of age, metallicity, internal extinction 
and mass can or cannot be disentangled in the various evolutionary stages throughout 
their lifetimes and what are the most useful passbands required to resolve the ambi- 
guities. We find the U and B bands to be of the highest significance, while the V band 
and near-infrared data provide additional constraints. A code is presented that makes 
use of luminosities of a star cluster system in all of the possibly available passbands, 
and tries to find ranges of allowed age-metallicity-extinction-mass combinations for 
individual members of star cluster systems. Numerous tests and examples are pre- 
sented. We show the importance of good photometric accuracies and of determining 
the cluster parameters independently without any prior assumptions. 

Key words: globular clusters: general - open clusters and associations: general - 
galaxies: star clusters - galaxies: evolution - methods: data analysis 



1 INTRODUCTION 

Since the seminal work by Tinsley (1968), evolutionary syn- 
thesis has become a powerful tool for the interpretation 
of integrated spectrophotometric observations of galaxies 
and galactic subcomponents, such as star clusters. Several 
groups introduced their evolutionary synthesis codes, e.g., 
Bruzual & Chariot (1993) [b&c], Leitherer et al. (1999) 
[starburst99] , Fioc & Rocca - Volmerange (1997) [pe- 
GASe], Fritze - v. Alvensleben & Gerhard (1994) [galev] (all 
with regular updates) , with various input physics (evolution- 
ary tracks vs. isochrones from various groups, different sets 
of stellar spectral libraries, extinction laws ...). The codes do 
not only vary in terms of input physics but also regarding 
computational implementation, interpolation routines etc. 
A number of publications deals with the intercomparison 
of various evolutionary synthesis codes (e.g. Worthey 1994, 
Chariot et al. 1996). The impact of uncertainties in the var- 
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ious model parameters (such as in the descriptions of over- 
shooting and mass loss, stellar spectral libraries etc.) on the 
resulting colours is challenged by Yi (2003). These publica- 
tions find a good general agreement among the various mod- 
els, and assign acceptable uncertainties to the model results. 
Yi (2003) points out the importance of a proper choice of 
filters for observing objects characterised by different age 
ranges. This is justified by the light being dominated by 
stars in different evolutionary stages at different times. The 
age-metallicity degeneracy is a major drawback for accurate 
age determinations, especially for young ages < 200 Myr. 

In addition to the choice of the specific evo- 
lutionary synthesis model used another important 
caveat merits discussion here. A common assumption 
in dealing with evolutionary synthesis is a well-populated 
stellar initial mass function (IMF), up to the model's upper 
mass limit. While this is probably a justifiable assumption 
for galaxy-sized systems (although uncertainties regarding 
the IMF slope persist), it certainly breaks down at levels 
of small (open) star clusters and OB associations, where 
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stars are formed purely stochastically (by consumption of 
the available amount of gas), and these statistics dominate 
the observed dispersion in cluster luminosities. A great deal 
of progress has been achieved already on this topic, in partic- 
ular by Cervino and collaborators (e.g. Cervino et al. 2002, 
Cervino & Valls-Gabaud 2003). The main conclusion is that 
for systems more massive than ~ 1O 5 M0 the impact of the 
stochasticity of the IMF on the results is - in general - low, 
and the UV continuum is least affected by stochastic disper- 
sions. 

The studies referred to before concentrated on 
the models themselves. When comparing the model 
results with observations, in order to constrain the 
cluster parameters age, metallicity, internal extinc- 
tion, and mass, one does not only need to take into 
account the model uncertainties, however. The final 
parameter uncertainties also depend on the obser- 
vational errors, the choice of passbands used, their 
number, spectral coverage and individual filter prop- 
erties, and the analysis algorithm applied to one's 
data. The most common way of model-observation 
comparison for astrophysical purposes is the chi- 
squared minimisation technique, used e.g. for pa- 
rameter determination of star clusters (e.g. Maoz 
et al. 2001, de Grijs et al. 2003a, b), determination 
of star formation histories of galaxies (e.g. Gavazzi 
et al. 2002), and photometric redshift determina- 
tion (e.g. Massarotti et al. 2001). Slightly different, 
but comparable algorithms, like the least-squares 
method (e.g. Ma et al. 2002) or maximum-likelihood 
estimation (e.g. Gil de Paz & Madore 2002, Bik et 
al. 2003), are used as well. However, see Bissantz 
&z Munk (2001) for a critical discussion about the 
applicability of chi-squared versus least-squares cri- 
teria. 

The aim of the present paper is a systematic 
evaluation of inherent uncertainties in the analysis 
of observed star cluster spectral energy distributions 
(SEDs) using evolutionary synthesis models. We de- 
fine an SED as an ensemble of (absolute) magni- 
tudes in a given set of (broad-band) passbands. We 
pay special attention to the most appropriate choice 
of passbands to improve future observation strate- 
gies. We will point out severe pitfalls, such as trends 
caused by finite observational errors and unjustified 
a priori assumptions. 



2 MODEL DESCRIPTION 

In section 12.11 we present the basic properties of 
our evolutionary synthesis models. Section 12.21 is a 
general description of our cluster SED analysis al- 
gorithm, regardless of whether it is used to study 
the parameters of observed star clusters or of simu- 
lated artificial clusters. In section !2.3l we present the 
specific properties of the artificial clusters (clusters 
for which SEDs are taken directly from our mod- 
els) used to simulate observed clusters and study 
the performance of our analysis tool. From section 
l3l onwards only these artificial clusters are used. 



2.1 Input Models 

We use the single stellar population (SSP) models presented 
in Schulz et al. (2002) , with important improvements regard- 
ing the treatment of gaseous emission in the early stages of 
the cluster evolution, as presented in Anders & Fritze - v. 
Alvensleben (2003). These models include isochrones from 
the Padova group including the TP-AGB phase, and model 
atmosphere spectra from Lejeune et al. (1997; 1998). These 
extend from 90 A through 160 /im for five different metal- 
licities, Z = 0.0004, 0.004, 0.008, 0.02 = Z Q and 0.05 or 
[Fe/H] = -1.7, -0.7, -0.4, and +0.4 (i.e., matching the 
metallicities of the Padova isochrones), and gaseous emis- 
sion (both lines and continuum) due to the ionising flux 
from young massive stars. The models can be retrieved 
from http://www.uni-sw.gwdg.de/~galev/panders/ For a 
general description of the stellar models see Bertelli et al. 
(1994) and Girardi et al. (2000), for details about the spe- 
cific isochrones in our models see Schulz et al. (2002). 

All calculations presented here are based on a Salpeter 
IMF in the mass range of 0.15 to approximately 70 Mq (0.15 
to approx. 50 Mq for super-solar metallicity; following from 
the Padova isochrones) . Stellar synthesis models for a Scalo 
IMF are presented in Schulz et al. (2002) and Anders & 
Fritze - v. Alvensleben (2003), and are available from the 
aforementioned WWW address. 



2.2 General description of the analysis algorithm 

In order to analyse observed SEDs of star clusters in terms of 
the individual cluster's age, metallicity, extinction, and mass 
we calculate a grid of models for a large range of values for 
each of these parameters (except mass, which is a sim- 
ple scaling of the model mass [M modcl = 1.6 x 1O 9 M ] 
to the absolute observed cluster magnitudes). Input 
parameters for the analysis are the time evolution of the 
spectra of the SSP models, and the derived magnitude evo- 
lution in the various passbands. 

The individual uncertainties contributing to the 
overall photometric uncertainties are: the observa- 
tional uncertainties, an estimated model uncertainty of 0.1 
mag, and an uncertainty of an additional 0.1 mag for pass- 
bands bluewards of the B band due to known calibration 
and model problems in the UV. The total uncertainty is the 
square-root of the quadratic sum of these individual errors. 
The observational and model uncertainties are expected to 
be independent. 

Galactic extinction is taken into account by dereddening 
the observations using the Galactic extinction values from 
Schlegel et al. (1998). 

First, we calculated dust-reddened spectra, using the 
starburst galaxy extinction law by Calzetti et al. (2000), 
assuming a foreground screen geometry, 

fe'(A) = 2.659 x (-1.857 + 1.040/A) + 4.05 
for 0.63 fj,m < A < 2.20/mi, 

fe'(A) = 2.659 x (-2.156 + 1.509/A - 0.198/A 2 + 0.011/A 3 ) + 4.05 
for 0.09 (im < A < 0.63/im 

with a reddened flux 

F Icd (\) = F (A) x 10 ^Vs(B-v)xk> W 
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and a range of values for the colour excess of the stellar 
continuum E s (B — V). Since the gaseous emission is relevant 
only for a short time and even then not the dominating 
term, the difference between the colour excess of the stellar 
continuum and that from nebular gas emission lines (e.g. 
Calzetti et al. 2000), is neglected. 

We emphasise that the Calzetti law is valid only for star- 
burst galaxies, while for "normal" galaxies (i.e., undisturbed 
and quiescent spiral and elliptical galaxies) it is probably at 
least marginally incorrect (due to the lower dust content 
in such galaxies). However, for our systematic uncertainty 
analysis, the specific shape of the extinction law assumed is 
of minor importance. 

We construct SEDs from these models by folding the 
spectra with a large number of filter response functions to 
obtain absolute magnitudes. The parameter resolutions 
are: 

• Age: 4 Myr resolution for ages from 4 Myr - 2.36 Gyr, 
20 Myr resolution for ages from 2.36 Gyr - 14 Gyr; 

• Extinction: the resolution is AE(_B — V) = 0.05 mag, 
for E(B -V) = 0.0 - 1.0 mag; 

• Metallicities: [Fe/H] = -1.7, -0.7, -0.4,0 and +0.4, as 
given by the Padova isochrones; 

• Mass: an arbitrary model mass of M mo del = 1.6 x 
1O 9 M is used. 

When comparing our observed SEDs with the 
model SEDs we first determine the mass of the clus- 
ter by shifting the model SED onto the observed 
SED. 

A number of these model SEDs (for Af c i ustGr = 
A/ m odci) are shown in Fig.[^ for the 5 available metal- 
licities and for 5 representative ages used for the ar- 
tificial clusters considered in this paper (see Section 

Each of the models in our grid is now assigned a cer- 
tain probability to be the most appropriate one, given by 
a likelihood estimator of the form p ~ exp(— \ 2 )i where 
X 2 — ( m ° bB ~y" modc1 ^ , where m b s and m m odci are the ob- 

°obs 

served and the model magnitudes in each band, respectively, 
and cr bs are the observational uncertainties. The summation 
is over all niters. Clusters with unusually large "best" x 2 are 
rejected, since this is an indication of calibration errors, fea- 
tures not included in the models (such as Wolf-Rayet star 
dominated spectra, objects younger than 4 Myr, etc.) or 
problems due to the limited resolution of the parameters. 
The cut-off level is set to a total probability < 10 -20 , corre- 
sponding to a Xbcst > 46. The total probability per cluster 
is then normalised. 

Subsequently, the model with the highest probability 
is chosen as the "best-fit model". Models with decreasing 
probabilities are summed up until reaching 68.26 per cent 
total probability (= 1 a confidence interval) to estimate the 
uncertainties in the best-fitting model. These uncertainties 
are in fact upper limits, since their determination does not 
take into account effects like the existence of several solution 
"islands" for one cluster (such as e.g. the age-metallicity de- 
generacy, see below), and discretisation in parameter space. 

For real observations, several passband combinations 
(containing at least 4 passbands) were used for the analysis, 
to minimise the impact of calibration errors and statistical 



effects. A minimum of 4 passbands is required to determine 
the 4 free parameters age, metallicity, extinction and mass 
independently (see also Anders et al. 2003, de Grijs et al. 
2003a,b). 

Only clusters with observational errors < 0.2 mag in all 
passbands of a particular combination are included to min- 
imise the uncertainties in the results (except for some arti- 
ficial clusters considered in this paper, for which we adopt 
errors = 0.3 mag). For each combination, the best-fitting 
models and their associated parameter uncertainties are de- 
termined. For a given cluster all best-fitting models (and 
the associated uncertainties) originating from the different 
passband combinations are compared. For each of these best- 
fitting models the product P of the relative uncertainties 

p agc^ ^ mass"^" ^ metallicity Z^" 

age - mass - metallicity Z — 

is calculated (the superscripts indicate the la upper ( + ) 
and lower (~) limits, respectively). The relative uncertainty 
in the extinction is not taken into account, since the lower 
extinction limit is often zero. The data set with the lowest 
value of this product is adopted as the most representative 
set of parameters (with its corresponding parameter uncer- 
tainties) for the particular cluster being analysed. In cases 
where the algorithm converges to a single model, a generic 
uncertainty of 30 per cent for all parameters is assumed, in 
linear space, corresponding to an uncertainty of \ 5 dex in 
logarithmic parameter space. See also Anders et al. (2003) 
for an application to the star clusters in the dwarf starburst 
galaxy NGC 1569, and de Grijs et al. (2003a,b) for applica- 
tions of this algorithm to clusters in the interacting starburst 
galaxies NGC 3310 and NGC 6745. 

2.3 Artificial clusters 

In this study we will use artificial clusters to investigate 
the uncertainties related to our analysis on the basis of a 
comparison with the model grid. The SED magnitudes 
of the "ideal" artificial clusters are taken directly 
from the models. Standard parameters of these clusters 
are: metallicity [Fe/H] = 0.0 = [Fe/H]©, internal extinction 
E(B - V) = 0.1, and ages of 8 Myr ("cluster 1"), 60 Myr 
("cluster 2"), 200 Myr ("cluster 3"), 1 Gyr ("cluster 4"), 
and 10 Gyr ("cluster 5"). In this standard set only a age 
variations, and neither metallicity nor extinction variations 
are considered initially, for reasons of clarity. The impact 
of varying the metallicity and extinction values is treated 
separately, see especially Sect. 13.31 The cluster mass is the 
model's mass 1.6 x 10 9 Mq, the "observational" errors are set 
to be 0.1 mag in each filter. Unless otherwise indicated, the 
clusters in this paper will have these standard parameters. 

For each of these 5 sets of artificial cluster parameters 
10,000 cluster SEDs were generated by adding statistical 
noise to the magnitudes of the "ideal" cluster. The errors 
are drawn from a Gaussian distribution with the Gaussian 
a corresponding to the "observational" uncertainty (= 0.1 
mag as standard value). 

All clusters are analysed separately with our algorithm 
in order to assess under which conditions and to what ac- 
curacy their input parameters are recovered by our method. 
Subsequently, all clusters originating from a given "ideal" 
cluster are used to calculate median parameters and their as- 
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sociated uncertainties. The uncertainties are centred around 
the median solution; they serve as equivalents to the la stan- 
dard deviation around the average values. However, for our 
analysis we chose to use the median instead of the average of 
the distribution, since we believe the median to be physically 
more relevant. We are interested in finding the most likely 
result when comparing our model grid with observations. 

Free parameters are the metallicity [Fe/H], the extinc- 
tion E(B — V), log(age) and log(mass). [Fe/H] and log(age) 
are used instead of Z and age because the evolution of mag- 
nitudes is approximately linear in [Fe/H] and log(age). 



3 STUDY OF THE ACCURACY OF OUR 
ANALYSIS 

3.1 Passbands included in our analysis 

We consider the following filters (the impact of only slightly 
different filter response curves is small) . All filters are taken 
from the set of available filters for observations of the Hub- 
ble Space Telescope (HSTJ/WFPC2, ACS, and NICMOS 
cameras. 

The standard set of filters is: HST WFPC2 (and ACS) 
filters F336W ("U"), F439W ( "B" ) , F555W ("V"), F675W 
("it"), F814W ("J"), NICMOS (NIC2 camera) F110W 
("J"), FI60W ("H"). This standard set will be referred to as 
"UBVRIJW . In addition the following filters are included 
in our study as well: the HST WFPC2 (and ACS where ap- 
propriate) wide filters F300W ("wide IF), F450W ("wide 
B"), F606W ("wide V"), F702W ("wide R") and the HST 
Stromgren filters F336W ("u" = "U"), F410M ("v"), F467M 
("b"), F547M ("y")- 

In this paper we will use the term "UV pass- 
band" essentially for the U band, and the term "NIR 
passbands" for the J and H bands. 

In the relevant figures, the horizontal lines mark the 
input values, and the symbols represent the median of the 
recovered values with the associated uncertainties. The clus- 
ters with "cluster number" = 1 < x < 2 are clusters with the 
youngest input age of 8 Myr, clusters with "cluster number" 
= 2 < x < 3 are clusters with an input age of 60 Myr, and 
so on (this offset is chosen for reasons of clarity) . 

3.2 Choice of passband combination 

First, we investigate which passbands contain the maximum 
amount of information, and hence which passbands are pre- 
ferred for observations, if one can obtain observations in only 
a limited number of passbands. This aims at improving fu- 
ture observing strategies. 

3.2.1 Importance of individual passbands 

In Fig. |5|we present the dispersions in our recovered param- 
eters using the standard input parameters, and SEDs cov- 
ering the full wavelength range UBVRIJH, compared with 
passband combinations where one of the UBVRIJH pass- 
bands is left out. 

This figure provides direct evidence of the importance 
of the U band (and to a lesser degree also of the B band) 
for all stages of cluster evolution, while for ages > 1 Gyr 



also a lack of the V band results in problems to recover the 
age. The systematic deviations from the input values for the 
combinations without the U or B bands are caused by an 
insufficiently accurate determination of the cluster metallic- 
ity. The resulting SED changes are therefore balanced by 
the analysis algorithm by adjusting the extinction and/or 
age, and are also accompanied by higher-than-input median 
masses in our fit results. 

Systematic biases are only apparent in the age deter- 
mination of the oldest artificial cluster (with a slight bias 
towards younger recovered ages), balanced by an overesti- 
mate of the internal extinction (which is a sign of the age- 
extinction degeneracy) and a minor bias towards smaller 
median masses. For the 60 Myr-old artificial cluster, the 
metallicity determination leads to an underestimate (pre- 
sumably due to the criss-crossing of the models and/or the 
non-negligible impact of the age-metallicity degeneracy at 
these ages) for all passband combinations, while for the old- 
est cluster the uncertainty in the metallicity determination 
encompasses almost the entire available range. 

In general, the median values recovered by our code 
agree fairly well with the input parameters, with the ex- 
ceptions mentioned above. The parameter dispersions are 
largest for the young (ages < 60 Myr) and the oldest (age = 
10 Gyr) clusters. This is caused by the criss-crossing of the 
models for young ages and the flat magnitude evolution for 
old ages. 

The importance of the U and B band is immediately ap- 
parent from the overview of artificial SEDs presented in Fig. 
Q U and B axe important for tracing the hook-like structure 
for young ages, while there appears to be a kink in the SEDs 
at the V band for older ages. 

3.2.2 Combinations of 4 passbands 

The minimum number of passbands required to determine 
the 4 free parameters age, metallicity, extinction and mass 
independently is four. 

In Figs. |3] and 0] we present the recovered parameters 
for UBVRIJH compared with various passband combina- 
tions consisting of 4 passbands, for optical filters only and 
including one near-infrared (NIR) band, respectively. 

For optical passbands only, the U band plays a ma- 
jor role once more, especially in determining the metallicity. 
Missing [/-band information leads to underestimates of the 
metallicity, thereby causing extinction values and ages to be 
adjusted improperly, and hence this also leads to incorrect 
mass estimates. Even in cases where the median is recov- 
ered correctly, such clusters show the largest uncertainties. 
In some cases, missing B-band information has similar ef- 
fects, especially for the youngest cluster, while for the old- 
est cluster the B band is vital to break the age-extinction 
degeneracy. Only for the oldest cluster, the V band contains 
vital information, which is in accordance with our results in 
Sect. 0231 

For optical+NIR passbands, the situation is similar: 
The U band (and to a lesser degree also the B band) is es- 
sential. Generally, the offsets from the input values and the 
uncertainty ranges are smaller than for optical passbands 
only, thus proving the importance of NIR data. Choosing 
a NIR band closely resembling the K band instead of J or 
H would give similar results, possibly restricting the values 
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Figure 2. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH and passband combinations rejecting 
one of the UBVRIJH passbands, as indicated in the legend. Cluster parameters are standard. 
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Figure 3. Dispersion of recovered properties of artificial clusters, assuming availability of various optical passband combinations, as 
indicated in the legend. Cluster parameters are standard. 
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Figure 4. Dispersion of recovered properties of artificial clusters, assuming availability of various optical+NIR passband combinations, 
as indicated in the legend. Cluster parameters are standard. 



slightly better. However, we concentrated on the H band 
since there are more observations available in H in the HST 
data archive than for filters with longer central wavelengths. 

In Fig.|3]we also see the effect of a limited wave- 
length coverage: in all parameters, the RIJH combi- 
nation gives the worst results (see also de Grijs et 
al. 2003a). Similar, but less pronounced is the effect 
for the UBVR combination. 

Fig.|S]compares the normal WFPC2 UBVR system with 
the corresponding passband combination using the WFPC2 
wide filters. In addition, results based on the medium band 
Stromgren filter system of WFPC2 are shown. 

In most cases the wide filter system gives slightly worse 
results than the standard system. However, driven by the 
wider filter response curves and the associated smaller ob- 
servational errors thanks to the larger flux throughput, the 
wide system might be preferable, e.g., for faint objects. 

Using the WFPC2 Stromgren medium-band system 
does not result in significant improvements compared to 
wide-band systems. In conjunction with the lower flux 
throughput (caused by the narrower bandwidth) this sys- 
tem seems less preferable for our purpose. We emphasise 
that this only holds for our SED analysis. 

In de Grijs et al. (2003a) we investigated the impact 
of the choice of passbands for the young cluster system 
(with ages of few x 10 - 100 Myr) in NGC 3310 with HST 
data from the UV through to the NIR. Starting with the 
full set of available passbands, we studied the changes in 
accuracy of the recovered parameters if we repeated the 
analysis using only a subset of our passbands. By compar- 
ing the results from our analyses using all passbands with 



those from smaller subsets we found severe biases in the age 
distributions originating from different passband combina- 
tions, in particular for combinations biased towards longer 
wavelengths (VIJH), but also for UV-UBV (covering shorter 
wavelengths only) and BVIJH, consistent with the results 
presented here. 



3.2.3 Conclusions on the choice of passbands 

From these comparisons we conclude that the pass- 
band combinations for the most reliable parameter 
determination must include the U band, the B band, 
and use the maximum available wavelength range, 
preferably including at least one NIR band. If only 
observations in 4 passbands can be obtained, the best combi- 
nations are UBIH or UBVH, especially for genuinely old ob- 
jects, and UBVI, if NIR data cannot be acquired. We empha- 
sise once again that tracing the kink around the B / V band 
in the SEDs (see Fig.^l is vital. For improved metallicity de- 
terminations, and consequently for improved determinations 
of the other parameters as well, NIR data seem to be cru- 
cial (for young clusters the U /B bands are also important, 
in order to determine the metallicity correctly). However, 
due to the limited metallicity resolution (and the numer- 
ous effects the metallicity has on the synthetic magnitudes), 
the metallicity determination remains the weakest point in 
our cluster analysis algorithm, and presumably in any rou- 
tine using synthetic magnitudes from stellar isochrones or 
tracks. 
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Figure 5. Dispersion of recovered properties of artificial clusters, comparing various wide and medium- band HST filters, as indicated in 
the legend. Cluster parameters are standard. 



3.3 Varying the input parameters 

In this section we investigate to which extent the input pa- 
rameters can be recovered as a function of their respective 
values and observational errors. 



3.3.1 Using all 7 filters 

Fig- 03 shows, for a range of observational uncertainties, the 
reliability of our recovered parameters if the standard set 
of filters (UBVRIJH) is available. We caution that we still 
apply the model uncertainty of 0.1 mag (and an additional 
uncertainty of 0.1 mag for UV passbands). 

A slight trend towards an underestimate of the ages, 
balanced by a slight overestimate of the internal extinction 
and an occasional underestimate of the metallicity, is seen. 
However, even for the largest observational errors of 0.3 mag 
that we tested for, all recovered parameters are consistent 
with the input parameters, within the uncertainties. 

With increasing observational errors, there seems to be 
a trend to underestimate the ages for the oldest cluster, bal- 
anced by an increasing overestimate of the internal extinc- 
tion. For genuinely old cluster systems, this degeneracy can 
be broken by restricting the extinction range. This is gen- 
erally justified, since such systems are usually dust-poor, if 
not dust-free, and show fairly homogeneous extinction dis- 
tributions. 

Fig.0shows that the degree to which our code recovers 
the input parameters is largely independent of the input 
extinction value, with the exception of the ages recovered for 
the oldest artificial clusters (in this latter case clear signs of 



the age-extinction degeneracy are apparent). The remaining 
deviations of the median recovered values from the input 
values are always less than 0.2 dex, and in most cases even 
smaller. The deviations in metallicity and extinction are one 
step in resolution (except for the extinction of the oldest 
cluster, which is, in most cases, 2 steps off). Small trends for 
increasing age underestimates with lower input extinction 
are discernible. 

Fig- HO indicates good agreement between the input pa- 
rameters and their recovered values for all 5 metallicities. 
Median extinction values and metallicities match the in- 
put values very well. The age determination is correct to 
A log(age) < 0.25 dex. The mass is recovered very well, as 
is the extinction. The various metallicity input values are in 
general correctly recovered, only in few cases a difference of 
one resolution step is seen. 

3.3.2 Using the minimum of 4 filters 

The following figures show the accuracy if observations in 
only the minimum of 4 passbands are available (i.e., a 
more realistic case). We discuss the best-suited 4-passband- 
combination identified in Sect. I3.2"!?l including the H band, 
i.e. the combination UBIH. 

Fig- HHshows significant trends caused by increasing ob- 
servational errors, especially for the oldest clusters. For the 
other clusters, the trends are less severe, with deviations of 
less than a factor of 2, or one step in the metallicity resolu- 
tion. For the oldest cluster, the deviations are up to 1 dex in 
age, 2 steps in metallicity, 0.35 mag in E(S — V) and 0.5 dex 
in mass, for the largest observational errors, i.e. 0.3 mag. 
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tional errors, as indicated in the legend. Other parameters are standard. 
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Figure 7. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH magnitudes and varying internal 
extinction values, as indicated in the legend. Other parameters are standard. 
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Figure 8. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH magnitudes and varying metallicity 
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Figure 9. Dispersion of recovered properties of artificial clusters, assuming availability of UBIH magnitudes and varying observational 
errors, as indicated in the legend. Other parameters are standard. 
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Fig. I1UI shows the recovered values for UBIH and var- 
ious input extinction values. For all but the oldest cluster, 
the recovered parameters reproduce the input values very 
well. The offsets and uncertainties are slightly larger, but 
comparable to the corresponding values for UBVRIJH. 

For the oldest cluster there are pronounced trends: with 
increasing input extinction, the recovered age, metallicity 
and mass estimates increase, while the offsets of the recov- 
ered extinction values from their input values decrease. We 
find that the higher the input extinction, the better all input 
parameters are recovered. 

Fig. llll shows the recovered values for UBIH and various 
input metallicities. The trends with increasing input metal- 
licity are less obvious, except again for the oldest cluster 
(more significant metallicity underestimates and extinction 
overestimates with increasing input metallicity). This be- 
haviour is also present, but less pronounced for the second 
oldest (= 1 Gyr old) cluster. For the other, younger clusters, 
the behaviour appears almost random. This is caused by the 
slow and steady evolution of magnitudes at large ages: with 
increasing metallicity the magnitudes become fainter and 
the colours redder. For younger ages the evolution of the 
magnitudes is less linear, and they criss-cross several times. 

3.3.3 Conclusions on the impact of varying the input 
parameters 

We have investigated the impact of varying the input pa- 
rameters on the accuracy of our parameter recovery. We find 
very good agreement, with generally small deviations, by ei- 
ther varying the input extinction or the input metallicity. 



Only for the oldest artificial clusters there are significant 
trends in the recovered values with increasing input extinc- 
tion and metallicity. This can be attributed to a number of 
degeneracies. We remind the reader that the upper 
age limit of the evolutionary synthesis models is 14 
Gyr. 

For increasing observational uncertainties there are 
clear trends of increasing recovered extinction and decreas- 
ing ages (with a small increase for the youngest cluster only), 
mass (an increase for the youngest cluster only) and metal- 
licity. The uncertainties increase as well, as expected. 

The results using either UBVRIJH or UBIH are fairly 
similar. Using 4 passbands only slightly increases the offsets 
of the median recovered values from the input values, as well 
as the uncertainties. Some trends, especially for the oldest 
cluster, become more significant. 



3.4 Restricting the parameter space to the 
correct ranges 

In this section we investigate the consequences of a priori 
restrictions of the parameter space. This might make sense 
in cases where, e.g., large observational errors may inhibit 
reasonable parameter constraints or where additional infor- 
mation is available, such as spectroscopic abundances, dy- 
namical age estimates for the starburst event that induced 
the cluster formation, a low-metallicity dwarf galaxy envi- 
ronment, etc. Here, we explore various cases where we re- 
strict some of our free parameters to the (correct) range 
of input values, and use passband combinations UBVRIJH 
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values, as indicated in the legend. Other parameters are standard. 
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and combinations lacking one of the UBVRIJH passbands 
to recover the input parameters. 

Fig. 1121 shows the results when we restrict the metal- 
licity to the (input) solar metallicity. Apart from the old- 
est cluster (which shows a slight underestimate of the age, 
balanced by a slight overestimate of the extinction) the re- 
covered values agree almost perfectly with the input values, 
in any case much better than without metallicity restriction 
(cf. Fig. 0. The importance of the 17 and B bands, and of 
the largest possible wavelength coverage are still reflected 
by the larger uncertainties for observations lacking those fil- 
ters. The deviations for the oldest cluster are a result of the 
age-extinction degeneracy. 

In Fig. 1131 we investigate the consequences of restrict- 
ing the analysis to the input extinction, allowing variations 
only in metallicity and age. The deviations of the median 
solutions from the input values, and the uncertainty ranges 
are significantly reduced compared to the unrestricted case 
shown in Fig. |2] Especially for the oldest cluster, some de- 
generacies are removed and the recovered values agree much 
better with the input values than in the unrestricted case. 
For genuinely old cluster systems the assumption of a generic 
low extinction may be justified, since such systems are com- 
mon in old relaxed galaxies with generally low (and uniform) 
dust content. 

Nevertheless, the importance of including the U and B 
bands is still apparent, especially in the age determination 
for the youngest cluster. By comparison with Fig. I12l we can 
attribute this behaviour to the age-metallicity degeneracy. 

The results from the restriction of both extinction and 



metallicity to their input values is presented in Fig. 1141 
Clearly, all median values agree perfectly with the values of 
the remaining input parameters. The few large uncertainty 
ranges indicate passband combinations that still do not al- 
low to determine the solutions unambiguously, because im- 
portant information (passbands) are missing. Combinations 
without these passbands do not allow for a reasonable anal- 
ysis. This includes combinations without the U or B band 
(and thus insufficient tracing of the kink in the SEDs) and 
without the U or H band (thereby restricting the wavelength 
coverage). The R/I bands seem to be of some importance in 
the early evolutionary stages to trace the curvature of the 
SEDs. 

Restricting the parameter space of our analysis to the 
input values for some parameters clearly reproduces the in- 
put values of the others, and hence can be used as a sanity 
check for the reliability of the algorithm. We find the age- 
extinction degeneracy to be most important for old clusters; 
for such systems a restriction in the allowed extinction range 
is usually possible. The age-metallicity degeneracy is respon- 
sible for some deviations for clusters younger than 200 Myr. 



3.5 Restricting the parameter space to incorrect 
values 

In this section we investigate the consequences of a priori 
assuming fixed, but incorrect generic values for the param- 
eters. 

First, we investigate the results for various input metal- 
licities, but using solar metallicity to recover the other input 
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Figure 14. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH and passband combinations 
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Figure 15. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH and various input mctallicities, 
as indicated in the legend. Solutions were sought with metallicity fixed to solar metallicity. Cluster parameters are standard. 
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Figure 16. Dispersion of recovered properties of artificial clusters, assuming availability of UBIH and various input metallicities, as 
indicated in the legend. Solutions were sought with metallicity fixed to solar metallicity. Cluster parameters are standard. 



parameters (as often done in the literature in studies of in- 
teracting and/or merging galaxies). The results are shown 
in Figs. ll5l and ll6l There are significant trends in the age de- 
termination in the sense of decreasing ages with decreasing 
input metallicity. These trends are in some cases accompa- 
nied by decreasing extinction. If the actual input metallicity 
is lower than the fixed metallicity assumed for the analysis, 
the cluster colours will be too blue for the combined input 
age and extinction, and for the fixed incorrect metallicity. 
Hence, either the recovered extinction is driven to lower val- 
ues and/or the solution to younger ages than their respective 
input values, since both adjustments result in bluer colours 
for ages > 200 Myr. The results for an actual super-solar 
input metallicity can be understood the other way around. 
The youngest clusters show rather randomly distributed re- 
covered values, thereby reflecting the complex behaviour of 
the magnitudes at such young ages. Applying solar metallic- 
ity models (as often seen in the literature) for clusters that 
intrinsically have sub-solar metallicity results in ages that 
are too low by up to 60 per cent, masses too low by up to 
56 per cent and similarly incorrect extinction values if ob- 
servations in 7 passbands are available, and even more if the 
observations were obtained in only 4 passbands. 

An equivalent test was done for our cluster sample in 
NGC 3310, where we compared the results from assuming 
a fixed, solar metallicity to leaving it as a free parameter 
(see de Grijs et al. 2003a). We found significantly different 
age distributions in either case, with the analysis in which 
the metallicity was left as a free parameter resulting in more 
realistic results in the context of what is known about the 
starburst in NGC 3310 in general. However, since all clusters 



are young in this cluster system (with ages of a few x 10 — 100 
Myr), this test was limited to young ages. 

In de Grijs et al. (2003b) we concentrated on the impact 
of restricting the allowed metallicity range for the analysis. 
By assuming a generic, fixed subsolar metallicity (confirmed 
by spectroscopy) , we found that the derived age distribution 
is fairly robust compared to the case where the metallicity is 
left as a free parameter, but the peak of the age distribution 
is significantly broadened. Hence, there is a larger dispersion 
for individual clusters. This is presumably caused by clusters 
that do not have the generic metallicity value. 

In Figs. [Inland CHI we show the results for UBVRIJH 
and UBIH, respectively, if clusters with various extinction 
values are analysed using a fixed value of E(_B — V). Shown 
are the results for clusters with input values E(B — V) = 
0.0, 0.1, 0.25, and 0.5, but analysed assuming a fixed extinc- 
tion E(_B — Infixed = 0.1. Considerable changes are observed 
in the resulting metallicities and ages. In many cases, the 
deviations from the input values are much larger than the 
derived uncertainties. For clusters with ages > 200 Myr 
there are significant trends of increasing recovered ages and 
metallicities with increasing input extinction. If the actual 
input extinction is smaller than the extinction fixed for the 
analysis, the cluster will be too blue for the combination of 
input age and metallicity, and for the fixed, incorrect ex- 
tinction. Hence, either the recovered metallicity is driven to 
lower values and/or the age must be younger than the cor- 
responding input values, since both result in bluer colours. 
The results for an actual input extinction higher than the 
fixed value can be understood the other way around. The 
youngest clusters show less obvious trends in the distributed 
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Figure 17. Dispersion of recovered properties of artificial clusters, assuming availability of UBVRIJH and various input extinction 
values, as indicated in the legend. Solutions were sought with extinction fixed to E(B-V)=0.1. Cluster parameters are standard. 
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Figure 18. Dispersion of recovered properties of artificial clusters, assuming availability of UBIH and various input extinction values, 
as indicated in the legend. Solutions were sought with extinction fixed to E(B-V)=0.1. Cluster parameters are standard. 
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recovered values with increasing input extinction, showing 
the complex magnitude behaviour for young ages. 

The results from this section again prove the im- 
portance of determining the physical parameters of 
the clusters, such as age, metallicity, and internal 
extinction (and mass as a derived value), indepen- 
dently, and avoiding any generic assumptions, which 
might not be justified for all clusters within a given star 
cluster system. This is true in particular for systems where 
the existence of two distinct cluster populations is already 
known or suspected, such as in merging galaxies and galaxies 
with known colour bimodality in their cluster system. 



4 CONCLUSION 

We have presented a detailed study of the reliability and 
limitations of our new algorithm to analyse observed SEDs 
of star clusters, based on broad-band imaging observations, 
by comparing these to a grid of model SEDs from our evo- 
lutionary synthesis code GALEV. 

We have computed a large grid of star cluster SEDs on 
the basis of our GALEV models for simple stellar populations 
including all relevant stellar evolutionary phases for ages 
> 4 Myr. The models also include metallicity-dependent 
gaseous line and continuum emission shown to be an im- 
portant contributor to broad-band fluxes in early evolution- 
ary stages. Our grid covers ranges in metallicity of —1.7 < 
[Fe/H] < +0.4, in extinction of < E(B - V) < 1, and ages 
of 4 Myr < age < 14 Gyr. The models produce spectra from 
which we derive absolute magnitudes, and hence broad-band 
SEDs, for any given filter system. Here, we present results 
for HST broad-band filters widely used for observations and 
analyses of star cluster systems in external galaxies. 

Our parameter analysis algorithm compares a given 
cluster SED (either observed or theoretical, as done in this 
study) with the model SEDs from our input parameter grid. 
Each parameter set is assigned a certain probability, based 
on an "observation-model" comparison using a chi-squared 
algorithm. The parameter set with the highest probability 
is adopted as the best model; the range of parameters from 
sets with the highest probabilities (up to a total probabil- 
ity of 68.26 per cent) determines the la uncertainties in the 
parameters. 

We constructed numerous artificial cluster SEDs, and 
varied each of the input parameters in turn to assess their 
effects on the robustness of our parameter recovery. For each 
clean model artificial cluster SED we calculated 10,000 ad- 
ditional clusters, with errors distributed around the input 
magnitudes in a Gaussian fashion. 

We identified useful and less suitable passband com- 
binations, with the aim to aid in the planning of observa- 
tional campaigns. Although a large number of passbands is 
always preferable, any realistic programme will more likely 
be limited to observations in the minimum number of re- 
quired passbands to successfully reach its goals. In order 
to successfully disentangle the three free parame- 
ters age, metallicity, and internal extinction based 
on the shape of a broad-band SED, and to deter- 
mine the mass of a star cluster by simple scaling of 
the model magnitudes to the observed level, a min- 
imum of four passbands are required. The most/least 



Table 1. Overview of the most important filters and most/least 
preferable 4-passband combinations, if NIR data are available 



Age 


important 


preferable 


combinations 




filters 


combinations 


to be avoided 


< few Gyr 


U, B 


UBIH, UBVH 


BVIH, RIJH 


> few Gyr 


B, V, U 


BVIH, UBVI 


UVIH, UBIH 



Table 2. Overview of the most important filters and most/least 
preferable 4-passband combinations, if no NIR data are available 



Age 


important 


preferable 


combinations 




filters 


combinations 


to be avoided 


< few Gyr 


U, B 


UBRI, UBVI 


BVRI, UVRI 


> few Gyr 


B, V, U 


UBVI 


UVRI, UBRI 



preferable passband combinations are summarised in Tables 
1 and 2 as a function of the expected age of the cluster 
population. In all cases, tracing the kink (or hook) 
in the SEDs around the B band (see Fig. [2> is of 
the highest importance. The inclusion of at least 
one NIR passband significantly improves the results, 
since NIR wavelengths allow to efficiently restrict the metal- 
licity range. For the youngest clusters, metallicity estimates 
are determined by the U and/or B bands. The poorest re- 
sults are obtained if neither UV information nor B band 
data are available, or if the available wavelength coverage is 
very short or biased towards blue or red wavelengths (like 
RIJH). 

By analysing artificial clusters, using a variety of input 
parameters (specifically age, metallicity, and internal extinc- 
tion) with our new code, we find in general good agreement 
between the recovered and the input parameters. Only the 
oldest, 10 Gyr-old artificial clusters show significant signs of 
the well-known age-metallicity-extinction degeneracy. 

We have considered several a priori restrictions of the 
parameter space, both to the (correct) input values and to 
some commonly assumed values. We easily recover all re- 
maining input values correctly if one of them is restricted, 
a priori to its correct input value; this also provides a san- 
ity check for the reliability of our code. We find the age- 
extinction degeneracy to be most important for old clusters. 
For such systems, an a priori restriction of the allowed ex- 
tinction range is often possible and shown to be very useful. 
The age- metallicity degeneracy is responsible for some mis- 
interpretations of clusters younger than 200 Myr. 

If we, however, restrict one or more of our input param- 
eters a priori to incorrect values (such as using, e.g., only 
solar metallicity, as often found in the literature), large un- 
certainties result in the remaining parameters. While certain 
restrictions might be justified in specific cases, we strongly 
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advice caution in more complex cases, such as in interacting 
galaxies or in galaxies with known colour bimodality in their 
cluster systems. 

Finally, we conclude that reliable determination 
of physical star cluster parameters is possible on the 
basis of broad-band imaging, provided the availabil- 
ity of a useful set of observational passbands, con- 
taining at least four filters, a sufficiently long wave- 
length base line, and reasonable photometric accu- 
racy. We show that a small, but suitably chosen filter 
set with deep observations (and the correspondingly 
small uncertainties) gives more reliable results than 
a larger number of shallow exposures in inappropri- 
ate or redundant filters. 

The method we have developed is a versatile and useful 
tool for the interpretation of large multi-colour data sets 
for star clusters of different ages and in a large variety of 
environments, such as provided by, e.g., our ST-ECF/ESO 
ASTROVIRTEL 1 project "The Evolution and Environmental 
Dependence of Star Cluster Luminosity Functions" (PI R. 
de Grijs). 
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